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ABSTRACT 

OO ' We investigate the recent and current star formation activity of galaxies as function of distance 

from the cluster center in a sample of 521 SDSS clusters at z < 0.1. We characterize the re- 
cent star formation history by the strength of the 4000A break and the strength of the Balmer 
\ absorption lines, and thus probe the star formation history over the last ~ 2Gyr. We show that 

when the Brightest Cluster Galaxies are excluded from the galaxy sample, there is no evidence 
for mass segregation in the clusters, so that differences in cluster and field populations can- 
<~i , not simply be attributed to different mass functions. We find a marked star formation-radius 

Qh' relation in that almost all galaxies in the cluster core are quiescent, i.e. have terminated star 

formation a few Gyr ago. This star formation-radius relation is most pronounced for low- 
Jl^ . mass galaxies and is very weak or absent beyond the virial radius. The typical star formation 

^ ' rate of non-quiescent galaxies declines by approximately a factor of two towards the cluster 

^ \ center However, the fraction of galaxies with young stellar populations indicating a recently 

completed starburst or a truncation of star formation does not vary significantly with radius. 
These results favor a scenario in which star formation is quenched slowly, on timescales sim- 
ilar to the cluster crossing time, i.e. a few Gyr. The fraction of star-forming galaxies which 

■ host a powerful optical AGN is also independent of clustercentric radius, indicating that the 
I link between star formation and AGN in these galaxies operates independent of environment. 
. The fraction of red galaxies which host a weak optical AGN decreases, however, towards the 

ff^ ■ cluster center, with a similar timescale as the decline of star forming galaxies. Our results 

can be fully explained by a gradual decline of star formation rate upon infall into the cluster, 
^k" . and rule out significant contributions from more violent processes at least beyond cluster radii 

^ ; > 0.17J200. 
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C3 ' 1 INTRODUCTION ellipticals are rarer in clusters jBamford et allllOOgl) . This suggests 



Galaxies in dense environments, i.e. in galaxy clusters and groups, 
are observed to have different properties than galaxies in less rich 



that star formation rate couples most strongly to environment, with 
morphology being only a secondary correlation. Evidence for this 
has been found in the local universe both in terms of local densit y 
environments - their star formation activity is subdued compared to jRauffmann et alJM Isianton et al.ll2005l:lBamford et aflB . 
field galaxies, and the occurrence of early-type galaxies is higher , ,■ r \ ■ " , ■ or,,,, A\^nnA 

I — --—I I 1 I ^ f r g ° as also rro m studies ot clusters (Christlein & Zabludonl 120051: 

than in the field dHubblel 1193a : IDresslen 119801 , and references |„ . , r , j t, u . i mno \-lt ■ \. 1 1 

lOuintero et al..,2006. : .van den Bosch et al..,200oa : .Weinmann et all 

thereto). In the field, star formation and morphology are not inde- ^QQgl^ 

pendent of each other: low star formation rate is correlated with 

early-type morphology (i.e. elliptical and SO galaxies), whereas tu- u-,- i» u , 

J \ f a jjj ^ hierarchical universe, clusters grow by accreting galax- 

high star formation activity is typical of late-type morphology (disk ■ r ,j ■ »Tj»r / ■■■■u 

° J jf f Bj \ les irom less dense environments. In order to torm / maintain the 

galaxies). This bimodality is modified in dense environments, as ^ r ^. j , .■ , .u j , u 

° ■' star formation - density relation, galaxy properties thus need to be 

evidenced by the atypical galaxies: the occurrence of passive spiral j ■ r ,, ■ , , » , r . . 

^ ■' *^ ° f r altered upon iniall into a group or cluster: star lormation needs to 

galaxies is higher in clusters than in the field, whereas star-forming , jj j-,, j.u^.j 

° ° "be suppressed, and some disk systems may need to be restructured 

to create spheroidal systems. Several mechanisms have been sug- 
gested by which either or both of th ese transitions may take place . 
* E-mail:anja@slac.stanford.edu There is increasing evidence (e.g. IWeinmann et al. I I2OO6I I2OO9I: 
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Ivan den Bosch et"aLll2008j |bl) that stranfjulati on, the stripping o f 
a diffuse gas reservoir surrounding the galaxy jLarson et al 1ll98(3t) . 
is the most important of these processes. Strangulation is effective 
already in the clu ster outskirts, and removes the gaseous envelope 
within a few Gyr l lBekki et al.ll2002h . With no diffuse gas available 
to replenish the cold disk gas, star formation is shut-off on sim- 
ilarly long time-scales. If the density of the intra-cluster medium 
is high enough, al so the cold disk gas may be stripped (ram- 
pressure stripping , lGunn&Gotllll972l : iFarouki & Shapiro! 1 1980l ; 
lOuilis et al.ll200ol) . thus truncating star formation on a very short 
time-scale. A few galaxies have been identified where ram-pressure 
stripping of the star-forming gas is evident through ta i ls of Hi or Ha 
emission and UV-bright kno ts I Gavazzi et al. 200 ll; Kennev et al] 



Hansen et alj|2009l) . as well as at higher redshift jPoggianti et al.l 
2006h . 



|2004 ISun & Vikhli"nQl2005l : ISun et al.ll2007|jcortese et alfeOOTl) 
all these galaxies are observed at the cores of massive clusters. At 
lower densities, i.e. in the cluster outskirts, ram-pressure stripping 
seems to be effective at truncating the cold gas disks as traced by 
HI beyond the stellar disk, but wi thout apparent effects on the opti- 
cal ap pearance of these galaxies l lCavatte et alj|l990l : lchung et al.l 



cat an 
l2007t) . 



By themselves, neither strangulation nor ram-pressure strip- 
ping can alter the morphology of infalling (spiral) galaxies, apart 
from the fading of star-forming regions and spiral arms. While 

galaxy mergers are an efficient mechanism to transform spiral 

' ' ' ^ ' 1 

galaxies into ellipticals (e.g. Toomre 1977; Farouki & Shapirg 

1982; Negrop on te & White 1983: Barn es & Hernquist 1992, 1996j; 
Naab & Burkertll2003l ; fCox et aflbooel) . they are unlikely to occur 



frequently in galaxy clusters due to the high relative velocities of 
cluster members (an exception is merging with the central galaxy). 
However, encounters between galaxies which do not lead to merg- 
ing can still affect galax y properties via tid al forces which heat the 
stellar disk (harassmen AMoore et alj[T998l) . Furthermore, galaxies 
experience tidal forces in the gravitation al potential of the cluster, 
which can similarly heat the stellar disk ( lGnediij|2003h . 

These processes can be distinguished by the different effects 
they have on the properties of galaxies, and by the timescale on 
which galaxies are affected. For example, strangulation causes an 
approximately exponential decline of the star formation rate, ram- 
pressure stripping truncates star formation abruptly, whereas merg- 
ers and tidal interaction cause gas to be transported to the center, 
thus triggering bursts of enhanced star formation. Examining the 
recent star formation histories of cluster galaxies is therefore a pow- 
erful tool to identify the processes at work. 

1.1 Radial cluster profiles as density probes 

In galaxy clusters, the star formation-density relation translates to 
a star formation-radius relation, since the density decreases with 
increasing (projected) clustercentric radius. Measurements of local 
density, such as the distance to the neighbor have the disad- 
vantage of strong shot noise, as well as selection and projection 
effects. In this sense, radial distance presents a much cleaner mea- 
surement. Stacking clusters allows for much better statistics than 
possible with any individual cluster, and yields the properties of 
an average cluster. Although the properties of individual clusters 
can differ significantly from the average cluster properties, even for 
clusters with similar masses (Moran et al. 2007), this scatter seems 
to be mostly stochastic in nature: several studies find that there is no 
systematic trend of varying gala xy populations with cluster mass, 
at leas t not for ma s sive clusters fenaka et al]|2004l ; iBalogh et al.l 
l2004bl ; lGotdl2005l ; iFinn et alj[200g) . There are, however, indica- 
tions of a dependence on cluster mass at lower masses (< IO'^^Mq, 



Clustercentric distance has a second interpretation apart from 
a local den sity measure: it also relates to the time since infall into 
the cluster jGao et alj2004l) . Of the galaxies at or beyond the virial 
radius, a considerable fraction has not (yet) experienced the dense 
cluster environment. Galaxies at the cluster core, on the other hand, 
were either born in a dense environment, or they must have tra- 
versed the cluster from the outskirts to the center at least once. The 
cluster crossing time is of the order of Rioq/cTv, i.e. 2.5 Gyr, and 
is independent of cluster mass, since i?200 cTv- Therefore, clus- 
tercentric distance is also an approximate timescale, and is sen- 
sitive to processes occurring on timescales of the order of Gyr. 
This timescale is particularly interesting for the processes which 
have been proposed to quench star formation in galaxies entering 
the cluster environment. Processes which quench star formation 
on long timescales, such as strangulation or harassment, would in- 
duce gradients over the complete radial range. Processes which act 
on shorter timescales are likely to cause distinctive signatures at 
the radii where they are most effective - e.g. ram-pressure strip- 
ping could be detectable as an increase in the (post-)starburst rate, 
presumably mostly in the cluster center, where the gas densities 
are highest. Mergers could induce (post-)starbursts primarily at 
outer cluster radii, where the relative galaxy velocities allow fre- 
quent mergers. (Note that the truncation of star formation induces 
the same spectral signatures as a true post-starburst galaxy, see 
Sect.[231) 



The dependence of the galaxy population mix on distance 
from the cluster center has been investigated in numerous stud- 
ies, with an overall consensus that the fraction of star-forming 
galaxies declines towards the cluster center, as traced both by 
galaxy co lor (e.g. 'Kodama & Bower'2001; Ouintero et al.l 120061 : 
iBlanton & B erlind 2007; Hansen et al. 20091) and emission-line 
strength (e.g. iLewis et alj|2 002: Go mez et alj|2003l : ISaiogh et al.l 
l2004al ; IChristlein & Zabludoff, .2005i) . These studies furthermore 
find that the suppression of star-formation can be traced to ~ 
2 - 4/?2oo- However, as we have motivated above, the properties 
of star-forming cluster galaxies can provide additional information 
that goes beyond a simple quantification of the galaxy bimodality 
in clusters. Studies of e.g. the occurrence of post-starburst galaxies 
have largely been li mited to intermedia t e-redshift clu sters, with the 
notable exception o f lHogg etal.l l [200^ . lHogg et al.l showed that in 
the local universe, post-starburst galaxies are fo und in the same 
enviro nments as star-forming galaxies (see also IZabludoff et al.l 
1 19961) . But they also find a marginal increase in the ratio of post- 
starburst to star-forming galaxies within the virial radius of clusters. 

At intermediate redshifts, the reported results are contradic- 
tory. Balogh et al. ( 1997, 199^) found no enhancement of starburst 
or post-starburst galaxies in a sample of 15 X-ray-selected clus- 
ters at z ~ 0.3 . But in a samp le of 10 clusters at s l ightly higher 
redshifts, [Dressier et akl i ll 9991) and IPoggianti etsi\ ( 1 199^ find a 
marked increase in the occurrence of post-starbursts. Subsequent 
studies have not been able to resolve these discrepant findings: 
evidence has been presented for clust er-triggered post-starbursts 
jTran et al.ll2004l ; |Poggianti et al.ll2009l) but also for suppres sionof 
post-starburst in groups jYan et aLlTiooi) . IPoggianti et alj ( |2009|) 
show that post-starbursts are found both in rich clusters as well as 
a in a subset of less massive groups. 
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1.2 The occurrence of AGN as function of environment 



1.3 This work 



A different viewing angle to the galaxy population mix is the oc- 
currence of Active Galactic Nuclei (AGN). The details of the link 
between AGN activity and the properties of the host galaxy, par- 
ticularly its star formation history, are still a matter of debate. But 
the mere existence of the linear correlation betwee n the mass of the 
central black hole and the mass of the bulge (e.g. iMagorrian et al.l 
1 19981) is a strong ind ication that the formation of the two mus t 
be closely hnked (e.g. lCattaneo et alj[T99^ ; lHeckman et al.ll2004l) . 
Powerful optical AGN trace current star formation activity, as they 
are fou nd predominantly in g alaxies with massive, but young, 
bulges jKauffmann et alj 1200 3l) . At least in the local universe, 
this AGN population correlates well with X-ray selected AGN 
jHeckman et al.l2005h . which are thus also likely linked to ongoing 
star formation. Weak optical AGN and radio-selected AGN, on the 
other hand, are fou nd predominantly in massive galaxies w ith old 
stellar populations jKauffmann et al.ll2003l : rBest et alj2005h . 

A prerequisite for AGN activity is a mechanism to transport 
gas to the central black hole. In disk galaxies, this may take place 
via major or minor mergers, or disk instabilities. Increased AGN 
activity in cluster galaxies could thus trace the effect of mergers, 
but possibly also harassment, which also calls for disturbance of 
the disk. 

However, the influence of environment on AGN activity is not 
yet well established. Large-scale statistical studies of optical AGN 
have become possible only with the SPSS, b ut a unanimous pic- 
ture has not emerg ed. iKauffmann et al.l ( 12004) showed that galax- 
ies in dense environments are less likely to host a powerful opti- 
cal AGN (L[Oni] > 10' Lq), similar to the suppression of star for- 
mation, but find no dependence on environment fo r weaker opti- 
cal AGN. A previous study by iMiller et alj l l2003h found no de- 
pendence of optical AGN occurrence on environment, presum- 
ably becaus e their sample was domina ted by weak AGN. On the 
other hand, IPopesso & Bivianol j2006l) report a lower fraction of 
(weak and strong) optical AGN in clusters with velocity dispersions 
greater than 500 km/s, compared to the field and smaller systems. 



In lBest et al.l l l2007h . we investigated the occuiTence of radio- 
selected and optical AGN in Brightest Cluster Galaxies (BCGs) 
compared to non-BCGs of the same stellar mass. We showed that 
BCGs are more likely to host a radio-loud AGN, but less likely 
to host a (powerful) optical AGN. We furthermore showed that 
enhanced radio-loudness occurs also for non-BCGs at the cluster 
core, whereas optical AGN emission is suppressed over most of the 
radial range to the virial radius. On the other hand, using a galaxy 

[ roup sample extending to much lower masses, IPasguali et al.l 
20091) find that both radio and optical AGN activity are enhanced 
in central galaxies vs. satellites of equal stellar mass, and that the 
occurrence of optical AGN is independent of clustercentric radius. 



In the local universe, X-ray selected AGNs are rare, and there 
is little evidence that the fraction of galaxies host ing an X-ray 
AGN differs systematically from clusters to the field l lMartini et al.l 
l2007h . There are however suggestions that the AGN fraction 
is hi gher in groups and p oor clusters compare d to richer clus- 
ters Icilmour et al] |2007| ; ISivakoff et all |2008|) . At higher red- 
shifts, there is evidence for a marked incre ase in the fraction of 
cluster galaxies wh ich host an X-ray AGN jEastman et al. I [20071 ; 
iMartini et af]|2009l) . but also that the AGN fraction is comparable 
between groups and the field (Georgakakis et al.„2008h . 



In this work, we investigate the composition of the cluster galaxy 
population as a function of clustercentric distance in a sample of 
521 clusters at z < 0. 1 from the SDSS. Sect.[2]presents the cluster 
and galaxy selection, as well as our diagnostics of the recent star 
formation history. We use spectroscopic indicators of current and 
past star formation rate, rather than just broad-band colors. We ap- 
ply new high signal-to-noise spectral indicators for light-weighted 
age and exc ess Balmer l ine str ength (to identify post-starbursts) de- 
veloped bv lWild et al. I l l2007l) . Using a diagram of light-weighted 
age vs. stellar mass, we divide the galaxies into stellar mass-limited 
subsamples of red (quiescent), green (transition-stage), young, and 
very young galaxies. In Sect.[3l we present the dependence of sev- 
eral galaxy properties, including nuclear activity, as function of 
clustercentric radius, and interpret these dependencies in light of 
the different galaxy transformation mechanisms. We summarize 
our results and conclusions in Sect.|4] 

Unless otherwise noted, we assume a concordance cosmology 
with Q.^ = 0.3, Ha = 0.7 and Ha = 100 /ikm/s/Mpc, where h = 
0.7. 



2 DATA 

2.1 Cluster sample 

Here we analyze 5 21 clusters from the sample defined in 

Ivon der Linden et al. I l l2007h . This cluster sa mple is derived from 
the C4 cluster catalog by Miller et al] | |2005|) . but special attention 
was paid in selecting the Brightest Cluster Galaxy (BCG) as the 
galaxy to be most likely to be at th e ce nter of the cluster potential 
well. In lvon der Linden et alj j2007h and lBest et al. 112003) we have 
shown that these galaxies indeed differ systematically from other 
galaxies. The differences can well be attributed to their special po- 
sition in the cluster center, and the rich merger history expected 
for central cluster galaxies. In the absence of X-ray centroids, the 
position of the BCG has been shown to be a good marker for the 
center, with a typica l offset from the X-ray centroid of < 100 kpc 
jKoester et"al]|2007h . We are therefore confident that the thus se- 
lected BCGs are the best available indicators of the cluster center 
for these systems. A further deviation from the original C4 catalog 
was the recalculation of the cluster velocity dispersions. 

For the present work, w e have further refined the sample from 
Ivon der Linden et al. I j2007h to account for the available data on 
cluster galaxies and the structure of the cluster For each cluster, we 
have det ermined the SDSS sur vey coverage from the photometric 
catalogs l lvon der Linden|[2007h . and retain only clusters for which 
more than 95% of the area within i?200 is in the survey area, i.e. 
is not subject to masks from bright stars, field edges, etc. Further- 
more, for each cluster we have determined a light-weighted cen- 
troid from galaxies in the photometric catalog, after doing a sta- 
tistical backg round subtraction to e stimate the likelihood of cluster 
membership ( I von der Linden|[2007l) . In the following, we consider 
only clusters for which the distance between the light centroid and 
the BCG is less than 400 kpc, and less than 0.3 R20Q ■ This con- 
straint is placed to ensure that the BCG position is a good tracer of 
the cluster center, and that therefore the association between dis- 
tance from the BCG and local density / time since cluster infall rate 
is valid. Clusters for which this is not the case, i.e. clusters with 
highly unusual structure, such as double clusters, are excluded from 
the analysis. These two constraints reduce the sample from 625 to 
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Figure 1. Distribution of the redshifts and velocity dispersions of the 521 
clusters in our sample. Error bars have been omitted for clarity. The his- 
tograms show the distribution binned in redshift (lower panel) and in ve- 
locity dispersion (right panel). The dotted lines indicate the divisons for the 
four redshift subsamples. 



585 and then to 521 clusters. The redshifts and velocity dispersions 
of these 521 clusters are shown in Fig.[T] 

2.2 Galaxy selection 

For the 521 clusters in our sample, we select the galaxies within 
±3cry and 20i?200i using the SDSS spectroscopic database. We ex- 
clude the BCGs from the galaxy sample, as we have shown that 
they differ systematica lly from other galaxies jvon der Linden et al.l 
l2007l : lBestetal]|2007h . The large radial range allows to probe ra- 
dial dependences to distances (far) beyond the turn-around radius. 
If a galaxy could be associated with more than one cluster, we as- 
sign it to the "closest" cluster, according to the distance measure 
A = '\j{d/R2oo)^ + (A,./(T,,p, where d is the radial distance from 
the BCG, and A,, the velocity offset. For this association of galax- 
ies with clusters, not only the final 524 clusters are taken into ac- 
count, but also those in the original C4 catalog. For those clus- 
ters that were removed in the BCG selection process, we use the 
velocity dispersion as given in the original catalog. According to 
iMiller et alT ( l2005h . the C4 cluster catalog is about 90% complete. 
Hence, our sample of field galaxies should be quite clean from 
galaxies in undetected clusters. Since the C4 catalog was compiled 
from DR3, only galaxies contained in DR3 are used for the field 
sample. Within i?200, these data are supplemented with DR4 data, 
for which the MPA/JHU value added catalogue^] were available at 
the time of writing. 

2.3 Spectral indices 

Our main focus in this analysis is the recent star formation history 
of cluster galaxies. Different spectral features are associated with 
different timescales of star formation. In this paper we make use 

' |http : //www . mpa-garching . mpg . de/SDSS/| 



of nebular emission lines which indicate star formation rates on 
timescales of 10' years, the Balmer absorption lines for timescales 
of ~ 0.5 Gyr and the 4000 A break strength for longer timescales of 
a few Gyr. For the lat ter two indicators, we make use of the spectral 
indices developed bv lWild et al. |(|2003), which are based on a Prin- 
cipal Component Analysis of the wavelength region 3750-4150A. 
The first principal component (PCI) is essentially equivalent to the 
strength of the 4000A break, and is thus a measure of the light- 
weighted mean age of the stellar population. The second princi- 
pal component (PC2) measures the excess Balmer absorption com- 
pared to the expectation value for a given 4000A break strength 
(the locus of galaxies with continuous star formation). Tradition- 
ally, the Balmer absorption line strength is quantified as the strength 
of the Kd line - the great advantage of the new index is that it can 
be applied to spec tra with consider ably lower SNR than required 
for measuring H(5 jWild et al.ll2007l) . 

Fig. m illustrates the distribution of a subsample of the SDSS 
galaxies used in this paper, in PC2 vs. PCI and the more traditional 
US vs. D„(4000) for comparison. On the right hand side of both fig- 
ures, the red sequence can be seen with large 4000A break strengths 
(PC1~ 0) from the predominantly old stellar populations. As we 
move left in both diagrams, the mean light-weighted stellar age de- 
creases, indicated by a decrease in 4000A break strength, and in 

this dataset a well defined blue sequence is seen at PCI 4. To the 

bottom left of both diagrams we find weak Balmer lines and weak 
break strengths i.e. starburst galaxies. And above the blue sequence, 
with strong Balmer absorption lines, we find the "Balmer strong" 
galaxies, which are often associated with post-starburst galaxies, 
or galaxies which have suffered a recent truncation in their star for- 
mation history. Because PCA identifies lines of correlation in the 
dataset, and 116 and D„(4000) are correlated in galaxy populations, 
the first component of PCA naturally takes account of this correla- 
tion. This causes the twist of the figures relative to each other, and 
means that PC2 is equivalent to "excess" Balmer absorption over 
that expected for the 4000A break strength of the galaxy. Along 
with the improved SNR of the PCA indicators, this makes the iden- 
tification of Balmer-strong galaxies easy, as they are displaced up- 
wards in PC2 from the main cloud of points. 

To aid interpretation of both diagrams we have overplotted 
examples of stellar population synthesis models (from Chariot & 
Bruzual 2009, in prep.), with differing star formation scenarios. The 
cyan points indicate the position of galaxies composed of the super- 
position of two components with differing star formation histories: 
an old "bulge" component, and a young "disk" component. Disk- 
to-bulge ratios vary between and 1, and the SFR of the disks is 
allowed to fluctuate on timescales of 10** years which causes the 
scatter in the distribution compared to standard smoothly exponen- 
tially declining SFHs (t models). The purple points are a subset 
of the cyan points - in these models, star formation has decreased 
substantially in a short timespan, i.e. has been truncated. Note how 
these star formation histories can easily be distinguished from the 
"normal" models by their larger PC2 values. The tracks show a dif- 
ferent star formation scenario: the superposition of short starbursts 
on old stellar populations. Galaxies then move in a clockwise circle, 
from the red sequence, quickly to the starburst phase, then climb in- 
creasingly slowly in PC2 or 116 into the post-starburst phase, until 
the burst fades and they reenter the red sequence. As we can see, the 
indices PC 1/PC2 are interesting for investigating the stellar popula- 
tions of cluster galaxies, since they can test whether cluster galaxies 
are prone to experience starbursts (e.g. through mergers) or trunca- 
tion (by ram-pressure stripping) more often than field galaxies. 
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Figure 2. In grayscale, the joint distributions of H5 vs. D„(4000) (left) and PC2 vs. PCI (right) for galaxies in our sample with masses < lO'"^ and redshifts 
> 0.06. The red sequence lies on the right of each diagram, and the blue sequence extends to lower D„(4000) or PCI. Starburst galaxies are found at the bottom 
left of each diagram, with weak Balmer absorption and weak 4000A break strengths. Post-starburst and truncated star formation galaxies lie above the blue 
sequence, with stronger Balmer absoiption than expected for their 4000A break strength. To aid in interpretation of the figure, we have created examples of 
population synthesis models with different star formation history scenarios. Overplotted as cyan dots are model galaxies composed of an old bulge population, 
and a young disk population with disk:bulge ratios varying from to 1 . The star formation history of the disks fluctuates around a constant, creating the scatter 
observed in the distribution. The purple points show models in which star formation has decreased substantially over a short timespan. This causes the models 
to move upwards in PC2 and diagonally to the top right in D„(4000)/H5. Overplotted as tracks are bulge+starburst models with burst mass fractions (defined 
as the total mass created in the burst divided by the total mass in the bulge) of 5% (red), 10% (blue) and 20% (green). Stars indicate time intervals 0.001, 0.01, 
0.1, 0.5, 1.0, and 1.5 Gyr after the start of the burst (the tracks run clockwise from the red sequence in both diagrams). 



2.4 PCI vs. stellar mass 

Fig. [3] presents the distribution of the galaxies in our sample 
in the PCI vs. stellar mass diagram in four redshift subsam- 
ples (see next section). The stellar masses have been determined 
from the model magnitud es, via the kcorrect_v4 .1.4 algorithm 
telanton & RoweislllOOTi) . Both a red sequence (strong 4000A 
break, PC1~ 0) and a blue sequence (PC1~ -4) are clearly visi- 
ble. 

This diagram is analogous to the more traditional color- 
magnitude diagrams (CMDs), in that both color and PCI are indi- 
cators of the recent star formation history, whereas luminosity and 
stellar mass measure the "size" of the system. PC 1 is a better indi- 
cator of the light-weighted age of the stellar population than broad- 
band colors, because it spans a much shorter range in wavelength, 
and is thus much less affected by dust extinction. Stellar mass is ar- 
guably a more fundamental galaxy property than luminosity, as the 
latter is highly dependent on the age of the stellar population: young 
stars are very luminous, hence galaxies with young stellar popula- 
tions have low MjL ratios. Furthermore, the MjL ratio varies with 
the observed band: blue bands are particularly dominated by young 
stars, whereas red bands are sensitive also to older stellar popula- 
tions. This means effectively that the red and blue sequences are 
shifted in CMDs of different bands. 

iBaldrv et al. I ( |2004|) fitted bimodal Gaussians to the red and 
blue sequences in a color-magnitude diagram of a large sample of 
SDSS galaxies in bins of luminosity. We follow the same approach 
here and fit the distribution in PCI in bins of stellar mass, shown 
in Fig. |4] Clearly, the distributions are not well described by two 
Gaussians - there is an excess of galaxies between the red and blue 
peaks. In CMDs, the region between the red and blue sequences 
has been coined the green valley, and we adopt this term here, as 



well. We fit two Gaussian curves to the distribution, but ignoring 
the green valley region for the fit. The fits provide adequate de- 
scriptions for the blue edge of the blue sequence and the red edge 
of the red sequence. However, in the green valley, there is clearly a 
population of galaxies not accounted for by these two components. 
The clear separation of the red and blue sequence, and identifica- 
tion of the green valley, illustrates the advantage of using spectro- 
scopic age indicato rs over broad-band colors, cf. Figs. 3 and 4 of 
lBaldrvetal.l ( l2004t) . 

Fig. |3] also indicates the peaks and Itr widths of the red and 
blue sequences. The ridgeline of the red sequence can be well de- 
scribed by two linear relations (the highest redshift bin does not 
probe low enough masses, hence only one line is fit), whereas 
the blue ridge is better described by a parabola. With these fits, 
we have continuous descriptions of the red and blue sequences 
with mass: PClRs(log M*/Mo) refers to the ridgeline of the red se- 
quence, PClBs(log M*/Mo) to that of the blue sequence. 



2.5 Galaxy classification 

The widths of both the red and blue sequences change with mass, 
most notably for the blue sequence (Fig.|4ll. This makes it unsatis- 
factory to define the extent of each sequence in terms of the widths; 
e.g. in the higher mass subsets, the green valley region is within 
1 - 2cr of the peak of the blue sequence, simply because the blue 
sequence is very broad. Instead, we use a fixed offset from each 
ridge to separate the three regions of the distribution. We define the 
upper (red) limit of the green valley to be half a unit in PCI lower 
than the peak of the red sequence, and the lower (blue) limit to be 
one unit larger than the peak of the blue sequence. As illustrated in 
Fig. [3] and Fig. |4] these limits correspond closely to the Icr width 
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z < 0.04 0.04 < z < 0.06 
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log(M./Ms) log(M./Ma) 

Figure 3. The distribution of PCI vs. stellar mass for all galaxies in our sample, in the four redshift subsamples, shown as logarithmically scaled grayscale. 
For each mass bin, the red (blue) circle shows the peak of the Gaussian curve fitted to the distribution of the red (blue) sequence, and is shown at the median 
mass of the particular mass range. The error bars indicate the respective Icr width. The dashed red and blue Unes indicate the Hmits of the green valley: galaxies 
above the dashed red line are classified as red, between the red and blue dashed lines as green, between the dashed and solid blue line as cyan, and below the 
solid blue line as blue (see Sect. 12. 5t . The vertical orange line indicates the mass limit for red galaxies in each sample. 



of both sequences for most of the mass range considered here. Fur- 
thermore, these limits reliably select the excess of green galaxies 
over the two gaussian fit in those mass bins where it is apparent 
(Fig.|4ll. By comparison, the Icr width of the blue sequence fails 
to do so at large stellar masses. Each galaxy is classified via the 
following scheme (also illustrated in Fig.|4l(: 

(i) If PCI > PCIrs - 0.5, the galaxy is classified as "red". These 
are the quiescent galaxies, in which star formation has ceased > 
1 Gyr ago. 

(ii) If PCIbs + 1 < PCI < PCIrs -0.5, the galaxy is classi- 
fied as "green", i.e. in the transition regime between young, star- 
forming and old, quiescent. 

(iii) If PCIbs < PCI < PCIbs + 1, the galaxy is classified as 
"cyan". Note that we have divided the young, blue population into 
two subsets - the "cyan" galaxies have slightly older populations 
than the "blue" galaxies. 

(iv) If PCI < PC1bs> the galaxy is classified as "blue". These 
are the youngest galaxies. 

While "green" may actually signify an intermediate class of galax- 
ies between the red and blue sequence, the term "cyan" has been 
introduced only to ease terminology, and "cyan" should be consid- 
ered a subclass of blue galaxies in the more general sense. Although 



PCI is strictly speaking a measure of the age of the stellar popula- 
tions, "blue" galaxies generally have higher levels of star formation 
than "cyan" galaxies. 

It is important to make the classification in this order: the red 
sequence is well defined in all mass bins (above the mass limit), but 
according to this definition, not every mass bin has a green valley. 
Extreme cases of this are the highest mass bins, where it is not 
possible to fit a peak to the blue sequence. 

The most important classification for the purpose of this paper 
is the distinction between old, quiescent (red) galaxies, and young, 
star-forming galaxies. We have verified that our results do not de- 
pend on the precise location of the cut between green and red galax- 
ies. In particular, all our results apply qualitatively also to the red 
side of the red sequence, i.e. galaxies redder than the red peak. 



2.6 Mass completeness 

Our aim is to work with a sample of galaxies complete in stellar 
mass, in order to follow the evolution of galaxies at fixed mass. 
The galaxies with the highest mass-to-light ratios are those with 
old stellar populations, i.e. the red galaxies . In Fig. |5]we plot the 
stellar masses of galaxies on the red sequence as a function of red- 
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mass range; 09.6-09.8 
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mass range; 09.8-10.0 



0.06 < z < 0.08 , mass range; 10.0-10.3 
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mass range; 10.5-10.7 
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mass range; 10.7-10.9 
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mass range: 10.9-11.1 
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mass range: 11.1-12.0 




Figure 4. The distribution in PCI for galaxies in the 0.06 < z ^ 0.08 cluster sample. The distributions are shown as black histograms. The green line shows the 
fit of a bimodal Gaussian to the data, where the region between the black dotted lines has been left out. The blue and red lines show the two individual Gaussian 
components. The short-dashed red and blue lines indicated the lo" widths of the red and blue sequence (only towards the green valley). The long-dashed red 
and blue lines show our adopted limits for the green valley - these are placed one PCI above the peak of the blue sequence and half a PCI unit below the peak 
of the red sequence. As can be seen e.g. from the 10.7 < log(M*/Mo) < 10.9 mass bin, the Icr width of the blue sequence falls into the PCI range with an 
excess of green galaxies, whereas the fixed-PCl-width cut correctly identifies this excess. The plot for the mass range 10.5 < logM© < 10.7 illustrates how 
these fits are used to classify the galaxies as blue, cyan, green, or red. Note that the mass ranges shown in the top two panels are below the mass completeness 
limit of red galaxies. 
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shift. For the whole sample, the sample is complete in red galaxies 
only for log(M<./Mo) > 10.5. But by defining subsamples of clus- 
ters (and their galaxies) in redshift bins, the completeness limit can 
be pushed to lower mass values. We define four such samples: the 
lowest redshift one (z < 0.04) is complete in red galaxies with stel- 
lar masses as low as log(M<-/Mo) > 9.6, but it contains only 39 
clusters. The sample with 0.04 < z < 0.06 (107 clusters) includes 
red galaxies with log(M,t/MG) > 10.0. The sample in the range 
0.06 <z< 0.08 (200 clusters) is complete for log(Mv,/Mo) > 10.3, 
and the highest redshift sample 0.08 < z < 0.10 (175 clusters) is 
complete only for log(M,,/Mo) > 10.5. 

Note that we refrain from using a V^ax weighting. This is be- 
cause Vmax assumes that galaxies are selected solely because of 
their apparent magnitude. This is not the case for our cluster galax- 
ies, as also the parent clusters are subject to selection biases. Fur- 
thermore, the effect of fiber collisions is not constant with redshift 
- at higher redshifts, fewer galaxies in denser environments can be 
targetted. 



cases where the light from the central 3" is not a good estimator 
of the total population. This would be the case e.g. for typical Sb 
galaxies, with on old, red bulge, and a young, blue disk. In the 
three redshift samples at z > 0.04, 98% of the galaxies identified 
as red from their PCI values also have a red overall color; in the 
z < 0.04 sample, this number is 90%. Furthermore, only a small 
fraction (~6%) of PCI -red galaxies have a dominant photometric 
disk component (frac_DeV< 0.5); of these galaxies, about 30-50% 
of these are indeed face-on, early-type spiral galaxies. The over- 
all "contamination" of the red sequence galaxies is therefore rather 
small, and of the order of 2-5%. We are thus confident that although 
the aperture bias distorts the PCl-mass diagram, our classification 
scheme is not significantly affected. 

The aperture bias further leads to an apparent shift of the loca- 
tion of the blue sequence with redshift, since the bias towards the 
central, oldest bulge population is largest in the closest galaxies. 
This is the reason why we need to fit the location of the blue and 
red sequence in separate redshift bins. 



2.7 Aperture bias 

In the PCI - stellar mass diagram, the red and the blue se- 
quences are not entirely parallel. At low stellar masses, the blue se- 
quence is approximately parallel to the red sequence. But at masses 
log(M*/MQ) > 10.3, the blue sequence curves upward towards the 
red sequence. This is a consequence of the aperture bias inherent 
to the SDSS survey, which obtains spectra within a 3" diameter. 
In our lowest redshift sample, 3" corresponds to ~ 2 - 3kpc, and 
thus only the inner regions of large (resp. massive) galaxies. Many 
massive disk galaxies also have a bulge component, with a stellar 
population older than the one in the disk. If only the inner part of 
the galaxy is observed, then the light of the bulge with its old stellar 
population is dominant. This is the reason the blue sequence turns 
"upward" with larger masses: generally, the bulge-to-total fraction 
increases with stellar mass for disk galaxies. Thus, the more mas- 
sive the galaxy, the more prominent the bulge, and the older the 
light-weighted population in the bulge. Note that in galaxy sam- 
ples based on slit spectroscopy, the red and blue se quences remain 
parall el and well-separated over all mass ranges l lvon der LindenI 
l2007h . 

The aperture bias could thus "contaminate" our sample in 



3 RESULTS: RADIAL PROFILES OF THE GALAXY 
POPULATION 

3.1 Mass segregation 

It has been claimed that the dependence of the galaxy popula- 
tion on environment is simply a consequenc e of a varying galaxy 
stellar mass function with l ocal density (e.g. iBamford et ah |2009| ; 
Ivan den Bosch et al .l2008d) . However, other studies based on clus- 
ter samp les find that any luminosity segregation is solely due to the 
BCG(s) l lBiviano et alj|2002l ; |Pracv et al.ll2005h . One would expect 
mass segregation in the clusters, if the dynamical friction timescale 
is comparable to the cluster crossing time. The more massive galax- 
ies then sink to the center faster than less massive galaxies, thus 
changing the mass function. In the left panel of Fig. [6] we test 
for this effect by e xamining the median mas s in the four redshift 
subsam ples. Unlike [Samford et al.l l i2009h and lvan den Bosch et al.l 
( l2008bl) . we find that the median mass appears constant with clus- 
tercentric radius. But the median as function of radius might be 
insensitive to changes in the radial distribution, and thus we test 
for mass segregation also via the cumulative radial distribution of 
four mass subsamples (right panel of Fig.[6ll. We do not combine 
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Figure 6. Left: The median mass as function of distance from tlie BCG in tlie four redsliift subsamples. To estimate the significance of each data point, 
we draw 100 bootstrap realizations of the cluster sample (rather than the galaxy sample). The error bars in this and the following figures indicate the 68% 
confidence intervals of the bootstrap. Right: The cumulative radial distributions within 2i?2oo of four mass bins in two redshift subsamples. Only the results for 
the two intermediate redshift samples are shown - the other two samples yield very similar results. Also note that the small offset seen between the two smaller 
mass bins and the two larger mass bins in the 0.04 < z < 0.06 redshift samples is significant only at the Icr level, and is not supported by the 0.06 < z < 0.08 
redshift bin, which has a much larger number of objects. 



the redshift subsamples in order to avoid any possible bias from 
fiber colhsions. If mass segregation were significant, it should show 
up as a deficit at small radii for the largest stellar masses. This is 
clearly not the case - the radial distributions for all mass subsam- 
ples are very similar. We thus find no indication of mass segrega- 
tion in our cluster samples. This likely indicates that the dynamical 
friction time in our clusters is larger than the characteristic cluster 
crossing tim escale for the majority of ou r galaxies. From N-body 
simulations. ISovlan-Kolchin et aLlfcoOSl) find that in order for the 
dynamical friction time-scale to be of the order of < 5 Gyr, the 
mass ratio between the satellite and the host halo needs to be large, 
> 0.1, and the merger orbit needs to be highly eccentric. For in- 
dividual galaxies and small groups infalling into the cluster, this 
is not the case, and the typical dynamical friction timescales are 
longer, > lOGyr The lack of noticeable mass segregation is thus 
fully consistent with such long dynamical friction timescales. 



How can this be reconciled with the results of Bamford et al.l 

( l2009h and lvan den Bosch et alj j2008d) ? lBamford et al.l do not ex- 
clude the BCGs from their study, and their s tellar mass trend ap- 
pears to be in fact largely driven by the BCGs. Ivan den Bosch et al.l 
do exclude the BCGs, but still find a significant trend of average 
stellar mass with radius. However, they find a strong gradient only 
for clusters of Migo < 10'''/7"'Mo. For larger halo masses, the gra- 
dient is significant only at the Itr level. Almost half of our clus- 
ters fall into this largest halo mass bin, and because these clus- 
ters are typically richer, they dominate our signal. The remaining 
discrepancy is thus not large. It is not clear what causes it - the 
cluster samples differ in many aspects. One possibility is that the 
discre pancy is caused by t he different adaptations of the cluster 
center: Ivan den Bosch et al.l measure the distance to the luminosity- 
weighted centroid, rather than the distance to the BCG. It is thus 



natural that they see luminous (and thus massive) galaxies at small 
radii. But particularly in case they missed the true BCG due to 
fiber collisions, or in the presence of significant substructure, the 
luminosity-weighted centroid may differ substantially from the po- 
sition of our BCG. 



3.2 Star formation-radius relation 

In Fig. [7] we investigate how each galaxy population (red, green, 
cyan, and blue) contributes to the number of galaxies within bins 
of mass and of clustercentric distance. For each mass range, we use 
the maximal number of clusters for which this mass range is com- 
plete, i.e. for the lowest mass range, 9.6 < log(Mt/Mo) < 10, only 
the lowest-redshift subsample with z < 0.04 is used, for the mass 
range 10 < log(M*/Mo) < 10.3 clusters with z < 0.06 are used (cf. 
Fig. |5)- For each mass range, we find a strong color-radius rela- 
tion: within about 1^200^ clusters are dominated by red galaxies. In 
the center, they make up more than 80% of the galaxy population, 
regardless of mass. The color-radius relation is particularly pro- 
nounced for low mass galaxies: in the lowest mass range, less than 
20% of the field galaxies are red galaxies. This fraction is approxi- 
mately constant till the virial radius, where it changes abruptly with 
a steep increase in the red fraction towards the cluster center. With 
increasing stellar mass, the red fraction in the field galaxies also 
increases, but also here, the strongest change is between ~ 0.3i?200 
and ~ 2i?2oo- 

By construction, the majority of field galaxies which are not 
red belong to the "blue" class. This is not the case within clusters. 
Clearly, the "blue" galaxies (i.e. the galaxies with youngest stellar 
population ages, and highest star formation rates) are most strongly 
affected by the cluster environment. The fraction of "cyan" galaxies 
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Figure 7. The fractions of red, "blue", 
"cyan", and "green" galaxies as function of 
distance from the BCG, in five mass ranges 
(indicated at the top of each panel). The 
color of the symbols and their connecting 
lines indicates the classification. 
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decreases more gradually. The fraction of "green" galaxies, how- 
ever, is almost constant across all radii. Only in the very cluster 
center (< 0.2/?20o) is it significantly l ower than in the field. T his is 
in good agreement with the results o f lWeinmann et al.l feOOd) . who 
found that the fraction of "intermediate" galaxies, as classified by 
g - r color, is independent of environment. 

These trends suggest that star formation is affected by a phys- 
ical process effective already at R200 and beyond. As star formation 
shuts off, "blue" and "cyan" galaxies redden and move onto the 
red sequence. This transition apparently does not occur instanta- 
neously, rather, these galaxies are seen as "green valley" objects for 
a time span somewhat shorter than a cluster crossing time. Hence 
the fraction of green galaxies decreases slower than the number of 
blue galaxies, as it is replenished by the latter ones as they redden. 
Similarly, "blue" galaxies will be first seen as "cyan", hence the 
fraction of the latter declines more slowly. 

One striking trend with mass is that the transition in the pop- 
ulation mix appears more abrupt at low masses, and more gradual 
at high masses. At low masses, suppression of star formation can 
traced to ~ 2i?20(h but at high masses, up to ~ 4/?2oo. We have veri- 
fied that this is not simply a consequence of cluster selection effects 
and the slightly different cluster samples employed for each mass 
range. One might conjecture that the more abrupt trend seen in the 
low-mass samples is simply a consequence of a possible higher 
completeness in selecting smaller galaxy groups at lower redshifts, 
and thus a purer field sample at low redshifts. The more gradual 
trend seen at higher stellar mass could then be caused by unde- 
tected galaxy groups in filaments around the cluster, which would 
already quench star formation prior to cluster infall. To test this 
hypothesis, we compared the star formation - radius relation in a 
fixed mass bin, 10.5 < logM© < 10.7, in all four redshift subsam- 



ples. The relation is remarkably similar in all four samples, and so 
we conclude that the qualitative differences in the mass subsamples 
is not caused by a possible redshift dependence of cluster selection. 

It has been suggested that an appreciable fraction of galax- 
ies obse rved at ~ 1 - 3i?200 previously pas s ed through the clus - 
ter (e.g. iBalogh et al.ll200(]| : lOill et alj|2005l : iLudlow et alj|2009h . 
This "backsplash" population should be dominated by low-mass, 
red galaxies. While we see suppression of star formation beyond 
R200, we see no evidence of this mass dependence. In all mass bins, 
the population at the virial radius and just beyond is very similar to 
the field. If anything, our data suggests that more massive galaxies 
are affected to larger distances from the cluster. 

The apparent influence of the cluster to distances well beyond 
R200 may be an indication of "pre-processing" of galaxies in the 
filaments and in small group s below our detection threshold (e.g. 
IZabludoflF & Mulchaevl[l998h . With increasing stellar mass, galax- 
ies cluster more strongly (e.g. iLi et al.ll200q) , and are thus more 
likely to be found in such environments. This provides a natural 
explanation for the earlier onset of the cluster influence for more 
massive galaxies. 

Finally, the sharp transition from field to cluster serves as ev- 
idence that our measured velocity dispersions and cluster centers, 
both of which required manual adjustments to an automatic pro- 
cess, are in general good representations of their true values. 



3.3 Star formation rate 

The decrease in the fraction of star-forming galaxies means that 
when averaging over the whole galaxy population, the star forma- 
tion rate clearly declines towards the cluster center. We have shown 
this via the ages of the stellar populations of cluster galaxies, which 
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Figure 8. The median star formation rate (left panel) and median specific stai' formation rate (right panel) of "blue", "cyan", and "green" galaxies which were 
classified as star-forming or low S/N star-forming galaxies. The dotted lines in the left panel show the median value for galaxies beyond the virial radius. 



probe the recent star formation (on timescales < 1 Gyr). If star for- 
mation indeed declines gradually in cluster galaxies, then also the 
instantaneous star formation rates in star-forming galaxies should 
decrease towards the cluster center. Thus, here we investigate the 
typical star formation rate of those galaxies with young stellar pop- 
ulation - the "blue", "cyan", and "green" galaxies. 

We use the sta r formation rate measurements of 
iBrinchmann et al ] (I2OO4I) . For consistency we use only galaxies 
which have been classified as star-forming or low S/N star-forming. 
The star formation rates for these galaxies are based on emission 
line model fits to the spectra. The star formation rates of other 
types of galaxies (AGN, composites, and unclassifiable) have been 
estimated from D(4000) and are thus representative of the recent 
star formation history (just as PCI) rather than the instantaneous 
rate. Note, however, that our results do not change if these galaxies 
are included. Just as the PCA, the star formation rate measurements 
are derived only from the central stellar populations covered by 
the fiber We show here the fiber star formation rates, but our 
results are qualitatively similar when using t he star formation 
rates c orrected for aperture effects provided by IBrinchmann et"al] 
( I2OO4I) . 

Fig. [8] shows the median star formation rate and the specific 
star formation rate (i.e. the star formation rate scaled by the galaxy 
stellar mass) of the star-forming galaxies as function of cluster dis- 
tance, again in five mass bins. The star formation rate clearly de- 
creases towards the cluster, though not very strongly (less than a 
factor of 2). This is consistent with the scenario we proposed ear- 
lier: as galaxies fall into the cluster, their star formation declines. 
The subsequent aging of the average stellar population causes the 
PC 1 value of the galaxies to increase - they transition from "blue" 
via "cyan" and "green" to "red". 

In each subpopulation of star-forming galaxies, however, the 
star formation rate varies very little with radius (Fig.|9j - it is pre- 
dominantly the proportions of strongly star-forming (i.e. young, 
"blue") galaxies, star-forming("cyan"), and weakly star-forming / 
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Figure 9. The median specific star formation rate of star-forming galaxies, 
for galaxies of stellar masses log(M^/Mo) > 9.6. ~ 7% of "red" ga laxies 
are classified as (low-SNR) star-forming bv lBrinchmann et al] )2004l) . 



transition ("green") galaxies which change and cause the decrease 
in overall star formation rate. 

The gradual decline of the typical star formation rate over 
most of the radial range from the cluster outskirts to the cluster 
center suggests that the quenching of star formation takes places 
on ti mescales similar to th e cluster crossing time, i.e. a few Gyr. In- 
deed, iBalogh et al.l ( |200(]|) modeled the star formation gradients in 
clusters, assuming that star formation declines on Gyr timescales, 
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Figure 10. The fraction of Balmer-strong "blue", "cyan", and "green" 
galaxies in two mass bins, as function of clustercentric distance. 



and predicted a decline by a factor of about 3, which is similar to 

what we observe. 

Our results differ from those o f lBaloghetal1 l l2004ah . who in- 
vestigated star formation rate as a function of environment density 
in SDSS and 2dFGRS galaxies. While they also found that the frac- 
tion of actively star-forming galaxies is a function of density, they 
did not find a significant difference in the distribution of Ha equiv- 
alent widths (as measure of specific star formation) of star-forming 
galaxies in different environments. A possible cause for this diflter- 
ence is that their density estimator, the distance to the fifth nearest 
neighbor, is a more noisy density measurement than clustercentric 
radius in our composite cluster. Furthermore, a pure local density 
estimator lacks a clear association with a time-scale since infall into 
a dense environment, which could be the determining factor. 



3.4 Fraction of Balmer-strong galaxies 

We have argued that the radial dependence of the central stellar 
populations of cluster galaxies suggests that star formation is not 
shut off instantaneously. This hypothesis can be tested via the oc- 
currence of galaxies with excess B aimer absorption line strength. 
If a galaxy experiences a brief starburst, or if a significant level of 
star formation is truncated on short timescales (< lO'yr), the galaxy 
spectrum displays strong Balmer absorption lines for ~ 0.5Gyr. 
The PCA measures the excess strength of the Balmer lines with 
the second Principal Component, PC2 (see Fig.O. Whi le the PCA 
is a p owerful tool to identify true post-starburst galaxies jWild et al] 
we are in this analysis interested simply in the extremes of 
the star-forming population. In each redshift and mass bin, we iden- 
tify the Balmer strong galaxies as the 5% of blue, cyan, or green 
galaxies with the largest PC2 values. This cut-off generally lies 
around PC2~ 0. Since PC2 is more sensitive to noise in the spectra 
than PCI, this analysis is carried out only for galaxies with SNR>8 
spectra (~90% of the galaxies). Because of the low number of ob- 
jects, we split the population only in two mass bins, with the divi- 
sion at the galaxy bimodality demarcation mass, 10'"^^ Mq. 

Fig. [To] shows the fraction of Balmer-strong blue, cyan, and 



green galaxies as function of distance from the cluster center. For 
low-mass galaxies, the fraction of Balmer-strong galaxies is con- 
sistent with being constant with radius. For high-mass galaxies, the 
data suggest a tantalizing increase in the two innermost bins, at 
< 0.2/?200- The signal is robust against variations in the definition of 
Balmer-strong galaxies, the signal-to-noise criterion, the mass bin- 
ning, etc. This is exactly the kind of signal one would expect from 
ram-pressure stripping, and it is tempting to associate it with this. 
However, the signal is carried by only 1 1 Balmer-strong galaxies in 
the first two bins. Furthermore, comparing the properties of the host 
clusters, these galaxies are found in clusters with similar velocity 
dispersions and X-ray luminosities as normal star-forming galax- 
ies, whereas t he efficiency of ram-p ressure should scale with cluster 
mass (but see lPoggianti et alj2009l for evidence for Balmer-strong 
galaxies in low-mass groups). Finally, the known galaxies in the 
process of being ram-pressure stripped tend to be low-mass galax- 
ies, but our data show no peak for low-mass galaxies. The increase 
in Balmer-strong high-mass galaxies at the core may therefore be a 
stati stical fluke. 

iHogg et"ai] ( l2006l) also found a marginal increase in the num- 
ber of post-starburst galaxies compared to the number of star- 
forming galaxies within clusters. They find an excess (compared 
to the field ratio) over all of the radial range up to R200, with the 
largest ratios at the cluster core (< 0.2/?2oo), as well as right at Riqq. 
While we cannot confirm the peak at Rim, nor an excess beyond 
0.2/^200, it is interesting that both studies find a peak excess at the 
cluster core, though it should be kept in mind that both are only 
marginal detections. 

With the possible exception of the cluster core, we find no 
strong environmental dependence of the ratio of Balmer-strong 
galaxies to star-forming galaxies. Since the fraction of star-forming 
galaxies decreases in the cluster, the absolute fraction of Balmer- 
strong galaxies also decrea ses towards the cluste r ce nter. This is 
consis tent with the results of jZabludoff et alj ( ll996h and lHogg et al.l 
(2006!), who found that post-starburst galaxies are predominantly 
field galaxies. In our terminology, this is simply a reflection of the 
higher field abundances of blue, cyan, and green galaxies. 

For low-mass galaxies, which are the dominant population of 
infalling, star-forming galaxies, the lack of a trend in the occurrence 
of post-starburst galaxies with cluster environment shows that the 
physical process(es) which cause a post-starburst signature in the 
spectrum operate no differently throughout most of the cluster en- 
vironment than in low-density environments. Since a (moderate) 
post-starburst signature results not only from a recent star-burst, 
but also from any rapid truncation of star formation, this indicates 
that the process(es) which quench star formation in these galaxies 
operate on longer timescales. 

Given that examples of (low-mass) galaxies in the process of 
being ram-pressure stripped have been observed, why do we not de- 
tect such a population statistically? One part of the puzzle is proba- 
bly that these objects are quite rare, as evidenced by the few objects 
known. One may speculate that it must take unusual circumstances 
for a small disk galaxy to travel to the cluster core with its star- 
forming disk intact. The other part of the puzzle is likely the selec- 
tion of clusters - our sample contains only few clusters of similar 
mass to the host clusters of these galaxies, and thus ram-pressure 
stripping is likely less efficient in our sample. 



3.5 Fraction of strong Active Galactic Nuclei 

iKauffmann eTal] l l2003h showed that powerful (L[Oiii] > 10^ L©) 
narrow-line AGN typically reside in massive galaxies with young 
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stellar populations, i.e. those galaxies which are the exception to the 
typical mass bimodality. Nuclear activity and star formation must 
therefore be tightly linked in these galaxies, and AGN could thus 
provide us with another channel to study environmental influences. 

A luminosity of L[Oni] > IO'Lq is approximately the demar- 
cation boundary between powerful Seyfert galaxies, and weaker 
LINERs. Furthermore, above this luminosity, AGN can be reliably 
identified in SP SS spectra regardless of redshift and star formation 
rate of the host jKauffmann et al.ll2003h . 

Fig. l3.5l shows the fractions of galaxies which host a powerful 
optical AGN with L[Ora] > IO'Lq. We see that the overall fraction 
of strong AGNs decreases towards the cluster center. Within each 
galaxy class, however, the fraction of strong AGN is consistent with 
being constant with radius, at least for radii > 0.2i?200- The decline 
of the total AGN fraction is then almost exclusively an effect of 
the increasing number of red galaxies, which host very few strong 
AGN. 

To increase the number statistics in the core of the cluster, 
we investigate the fraction of strong AGN in blue, cyan, and green 
galaxies combined (the purple symbols in Fig. l3.5l l. We do not find 
significant evidence for an enhanced fraction of strong AGN even 
in th e cluster core. 

iLi et alj ( l2008h and lReichard et alj ( |2009|) showed that kine- 
matical disturbances such as close companions, mergers, and disk 
lopsidedness trigger star formation at the centers of galaxies. While 
they also find an increased occuiTence of optical AGN, this is en- 
tirely acco unted for by the existence of th e young stellar popula- 
tions. Also lKauffmann & Heckmai] ( |2009|) argue that for the pow- 
erful AGN hosted by star-forming galaxies, black hole accretion 
rate is largely independent of processes in the interstellar medium, 
as long as the central stellar population is young. These studies 
suggest that the star formation rate at the galaxy center regulates 
the probability that the central black hole shines as a powerful op- 
tical AGN. Our results are consistent with this scenario: our galaxy 
classification scheme is based on the age of the central stellar pop- 
ulation, and thus the AGN fraction within each galaxy class is inde- 
pendent of cluster radius. Due to the suppression of star formation, 
the overall AGN population is suppressed in the cluster environ- 
ment, as well. Strictly speaking, this scenario precludes the use of 
powerful AGN as tracers of environmental processes. 

3.6 Fraction of AGN in red galaxies 

While the more powerful Seyferts are found almost exclusively in 
star-forming galaxies, the weaker LINERs are found predominantl y 
in galaxies with old stellar populations jKauff^mann et al.ll2003h . 
The star formation rate in our red galaxies is quite low (Fig. [gj, 
so that AGN can be identified to much lower luminosities without 
contamination from star foiTnation emission lines. Fig. 1121 shows 
that the fraction of optical AGN in red galaxies declines towards 
the cluster center, by approximately a factor of two between the 
field and the center. The decline is gradual, and consistent with set- 
ting in around Rjoo- 

The relation between these weak AGN and their host galaxy 
properties is very different fro m the more powerful optical AGN. 
iKaufTmann & HeckmarJ | |2009|) showed that for these systems, the 
rate of accretion onto the black hole is proportional to the bulge 
mass, with a constant of proportionality that depends on the age 
of the bulge stellar population. A likely explanation for this re- 
sult is that the AGN is fuelled by mass loss from evolved stars. 
D„(4000) and PCI are largely insensitive to stellar populations of 
> 2 Gyr, and thus our class of red galaxies is far from a homo- 
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Figure 12. The fraction of red galaxies which host an optical AGN. 



ge neous population. It is ther efore plausible (and indeed expected, 
see lPe Lucia & Blaizoll2007l) that the mean age of the red galaxies 
increases towards the cluster center, and provides less fuel for the 
central black hole. 

There is furthermore evidence that the (weak) AGN activity in 
galaxies with old central stellar po pulations may be sensiti ve to pro- 
cesses in the outer galaxy regions. iKauffmann et al.l ( l2007h find that 
galaxies with an old central stellar population which host an optical 
AGN span a wide range in UV-optical color. In contrast, galaxies 
with old central populations, but no AGN emission, almost all have 
red UV-optical colors. They interpret the UV emission as residual 
small levels of star formation in an outer disk, traceable only via ul- 
traviolet radiation. These star-forming outer disks presumably trace 
a (low-mass) reservoir of cold gas, which is available to fuel star 
formation and AGN activity^ Given that we find that fewer red 
cluster galaxies host an AGN, we expect that the fraction of red 
galaxies with a UV-bright outer disk also decreases in the cluster. 
Just as star formation in young galaxies is suppressed upon infall 
into the cluster environment, also this "hidden", residual star for- 
mation in otherwise old galaxies is quenched, presumably because 
potential fuel is removed from the outer regions. 



4 SUMMARY AND DISCUSSION 

Our results can be summarized as follows: 

• The median mass of cluster galaxies other than the BCG 
does not vary with cluster radius. Environmental differences in the 
galaxy population are therefore not solely a result of a varying stel- 
lar mass function. 

• There is a pronounced relation between distance from the clus- 
ter center and the composition of the galaxy population. In the clus- 
ter center, star formation has ceased in ~80% of the cluster galax- 
ies, regardless of stellar mass. The fraction of the youngest ("blue") 

^ Note, however, that this is not a linear relation. Galaxies with UV-blue 
outer disks span a wide range in central star formation rate, age of the cen- 
tral stellar population, and black hole accretion rate. 
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Figure 11. The fraction of powerful optical 
AGN (L[Oiii] > IO'Lq) as function of clus- 
tercentric distance, in each of the five mass 
bins. The black dots and lines show the over- 
all fraction, the red, green, cyan, and blue 
Unes show the fraction of the con'espond- 
ing galaxy populations. The purple symbols 
show the fraction in "green", "cyan", and 
"blue" galaxies combined. 
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Stellar populations declines most rapidly towards the cluster center, 
less rapidly for slightly older ("cyan") stellar populations, and is 
almost constant for the transition "green valley" galaxies. 

• The star formation rates of galaxies which are not on the red 
sequence declines by about a factor of two from the field to the 
cluster center. 

• The fraction of "blue", "cyan", and "green" galaxies which 
have excess B aimer absorption line strength, indicative of a post- 
starburst or truncation of star formation, is independent of cluster 
radius at least for > 0.2i?200- There is a marginally significant in- 
crease of Balmer-strong galaxies within 0.2/?2oo for galaxy masses 
> 10'°-^ Mg. 

• The fraction of star-forming and transition galaxies which host 
a powerful optical AGN (L[Oiii] > 10' Lq) is independent of cluster 
radius. 

• The fraction of quiescent (red) galaxies which host a (weak) 
optical AGN decreases towards the cluster center. 

Our results indicate that the star formation rate of galaxies 
infalling into clusters decreases gradually over a timescale simi- 
lar to the cluster crossing time, i.e. a few Gyr. E.g. a simple ap- 
proximately exponential decline of the star formation rate suffices 
to explain the declining fraction of very young (i.e. strongly star- 
forming) in the cluster. These "blue" galaxies would then be seen as 
star-forming, "cyan" galaxies and later as transition,"green" galax- 
ies, which naturally explains the slower decline of these two types. 

Remarkably, the properties of star-forming and transition 
galaxies, at fixed stellar mass and fixed light-weighted age, are 
largely independent of clustercentric radius. The star formation rate 
in each of these classes is nearly constant, as well as the frac- 
tion of galaxies which host a post-starburst galaxy or a strong op- 
tical AGN. This argues against more violent processes than the 



mere fading of star formation being a dominant effect: truncation 
of star formation due to e.g. ram-pressure stripping would lead to 
an increase of Balmer-strong galaxies, whereas mergers and kine- 
matic instabilities (caused e.g. by harassment) would lead to an in- 
crease of starburst galaxies and Balmer-strong galaxies, and pos- 
sibly AGN activity. There is no evidence of any of these over the 
radial range we can probe, apart from a marginally significant ex- 
cess of Balmer-strong galaxies the cluster core. 

We therefore come to the conclusion that a simple gradual de- 
cline of star formation suffices to explain our results, and that other, 
more abrupt, processes cannot play a major role at cluster radii be- 
yond a tenth of the virial radius, at least not in the moderately mas- 
sive clusters we probe here. 

Our data alone can only approximate the timescale on which 
star formation is quenched in clusters - it must be of the order of 
the cluster crossing time, a few Gyr In Weinmann et al., in prep, we 
will compare a similar dataset to semi-analytical models in order to 
investigate the timescale more precisely. 
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